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Abstract 

Wc investigate the stability of black hole accretion disks in a primordial environment (POP III disks 
for short), by solving the vertical structure of optically thick disks, including convective energy transport, 
and by employing a one-zone model for optically thin isothermal disks. Because of the absence of metals in 
POP III disks, wc find significant differences in stability associated with ionization between POP III disks 
and the disks of solar metallicity. An unstable branch in S-shaped equilibrium curves on the 7\/ — E (mass 
accretion rate - surface density) plane extends to a larger surface density compared with the case of disks 
of solar metallicity. The resulting equilibrium loci indicate that quasi-periodic oscillations in luminosity 
can also be driven in POP III disks, and their maximal luminosity is typically by an order of magnitude 
larger than that of the disks of solar metallicity. Such a strong outburst of POP III disks can be observed 
by future huge telescopes, in case that the mass is supplied onto the disks at the Bondi accretion rates in 
typical virialized small dark halos. 
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1. Introduction 

According to recent advances in the formation theory 
of population III (POP III) stars, it has been understood 
that they could be very massive (10^ — IQ^Mq) (Omukai 
& Nishi 1998; Abel et al. 2000; Bromm et al. 2002). 

Those stars end up as black holes as massive as 
10- IO^Mq (here after POP III BHs), in the case that 
the mass of the progenitor stars is within the range 
40A/o<M* <140Afo or Mh.>260Mo (Heger & Woosly 
2002). The accretion disks formed around such black 
holes can emit a significant amount of soft X-rays, which 
can contribute to reionization of the universe (Ricotti & 
Ostriker 2004). Moreover, such black holes could be candi- 
dates for the progenitors of ultra-luminous X-ray sources 
(Mil & Totani 2005) or the building blocks of super- 
massive black holes (Rees 1984). 

The accretion disks surrounding POP III BHs are 
formed from the accreting gas around the progenitor stars. 
In the case that the mass of a POP III star is within the 
range AOMqKM^KIAQMq or >26OM0, the star forms 
a black hole directly without a SN (Heger & Woosly 2002). 
In this case, the accreting gas is expected to contain a 
small amount of heavy elements. In addition, recent nu- 
merical simulations on low-mass galaxy formation, taking 
into account the effects of SNe feedback, suggest that the 
metallicity of the gas remain in such galaxies during their 
early star-burst phase is as low as ~ 10~^Zq (Wada & 
Venkatesan 2003). Therefore, the metallicity of the ac- 
cretion disks surrounding POP HI BHs is expected to be 
very low {Z <1Q~'^Zq). We hereafter call these accretion 
disks as POP HI disks. 



The thermal instability of the standard accretion disk 
has been investigated in detail (e.g, Shibazaki & Hoshi 
1975; Shakura & Sunyaev 1976; Pringle 1976). By assum- 
ing solar metallicity (hereafter POP I disks), Mineshige 
and Osaki (1983) obtained the S-shaped equilibrium loci 
on the M — S (mass accretion rate v.s. surface density) 
plane, implying that the unstable disk undergoes limit- 
cycle oscillation between the upper, hot, ionised branch 
and the lower, cool neutral branch (see also Hoshi 1979; 
Meyer & Meyer-Hofmeister 1981; Smak 1982; Kato et al. 
1998 for a review). It is so-called dwarf- nova type instabil- 
ity, which successfully explains the observed properties of 
dwarf nova (Meyer & Meyer-Hofmeister 1984). An X-ray 
nova is also thought to be caused by such a disk instabil- 
ity, although its central object is expected to be a black 
hole (Huang & Wheeler 1989; Mineshige & Wheeler 1989). 
The limit-cycle behavior for active galactic nuclei has been 
investigated by Lin & Shields (1986) (see also Mineshige 
& Shields 1990). 

Mayer and Duschl (2005) investigated this instability on 
the POP HI disks in which they also obtained the equilib- 
rium loci on the M — S plane. However, their results based 
upon a one-zone approximation, even if the disk is opti- 
cally thick, in which the energy transfer due to convection 
along the direction perpendicular to the equatorial plane 
of the disk is not included. As a result, the equilibrium 
loci are not S-shaped, which does not directly lead to the 
limit-cycle oscillations driven by thermal and secular in- 
stabilities. Thus, we need a more sophisticated approach 
in order to understand the nature of the stability of POP 
HI disks. In this paper, we re-examine stability of POP III 
disks, by solving the vertical structure of optically thick 
disks, including the effect of convection. For optically thin 
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disks, we adopt a one-zone approximation. In addition, wc 
evaluate the amplification of the luminosity of accretion 
disks surrounding POP III BHs by the limit-cycle oscil- 
lation. Wc also discuss the possibility to detect the POP 
III disks by calculating their emergent spectra in the burst 
phase. 

This paper is organised as follows: we describe basic 
equations in section 2. We compare the stability of POP 
III disks with POP I disks in section 3. In section 4, we 
discuss the observational possibility of POP III black hole 
accretion disks, and summarise in the final section. 

2. Model and Numerical Method 

2.1. Optically Thick Regime 

In this subsection, we show the basic equations used to 
describe the vertical structure of the optically-thick ac- 
cretion disks. Here, we assume that the disks are in hy- 
drostatic equilibrium along the vertical axis, and local en- 
ergy balance is attained. The convcctive energy transport 
is taken into consideration. The equation of hydrostatic 
equilibrium is given by 



dp 
dz 



GMz 



(1) 



where p is the total pressure, p is the density, M is the 
black hole mass, and r and z are the radial and vertical 
coordinates. The equation of energy transport is 



dhiT 
d\nz 



Vrad 
^conv 



(V,ad < Vad) 
(V,ad > Vad) 



(2) 



where T is the temperature, and Vad is the adiabatic gra- 
dient. 

The radiative gradient, Vradi is written as 



Vrad — 



'Sr^KpF 



(3) 



AacT^GMz ' 

where n represents the opacity, F and a denote the en- 
ergy flux and the radiation constant, respectively. We 
utilise the latest OPAL opacity tablc^ including low tem- 
perature opacities (Ferguson et al. 2005). The energy flux 
is evaluated by integrating the following equation, which 
represents the local energy balance : 



dF 3 

1^ = -r^"- 



(4) 



Here a is the viscosity parameter, whereas — 
(GM jr^Y^"^ represents the Keplerian angular speed. We 
employ the expression ol = min[1.0, 10^(iJ/r)^-^] (with H 
being the disk half thickness), which succeeds to reproduce 
the variation amplitude and the timescale of the X-ray no- 
vae (Kato et al. 1998, see also Meyer & Mcycr-Hofmeister 
1984), whereas the viscosity model is still poorly under- 
stood for the primordial composition. We here note that 
the source of disk viscosity is though to be of magnetic 
origin based on the recent magnetohydrodynamic simu- 
lations (Hawley et al. 2001; Machida et al. 2001; Balbus 
2003, for a review). 



The convcctive gradient, Vconv, is calculated by the 
mixing-length formalism for convection (Paczyhski 1969), 
in which we assume mixing length, to be 



; = min(i?,i/p). 



(5) 



where Hp is the pressure scale height. 

In order to close the set of equations (l)-(5), we need the 
equation of state for the system in radiative equilibrium, 



pkT 1 



(6) 



/imp 3 

where nip denotes the proton mass, k represents the 
Boltzmann constant, and /i is the mean molecular weight. 

We integrate equations (1), (2), and (4) from the disk 
surface {z = H) toward the equatorial plane (z = 0) by the 
Runge-Kutta method for given parameters, r, Z, and M, 
where Z is the metallicity and M is the mass accretion 
rate. The outer boundary (z = H) conditions for temper- 
ature, energy flux, pressure are given by 

-(^^)"". 

F^aT^, (8) 



and 



P- 



3k 



(9) 



where a denotes the Stefan-Boltzmann constant. Wc it- 
eratively search H so as to meet the condition F = at 
z = 0. Then, we obtain the vertical structure of the ac- 
cretion disks. In this method, the column density and the 
optical thickness are calculated from E = 2 pdz and 

r = pndz, respectively. This method breaks down in 
the regions of a small mass accretion rate for POP HI 
disks, since then the disks become optically thin (We will 
discuss this point later) . A numerical method for optically 
thin disks is described in the next subsection. 

2.2. Optically Thin Regime 

The optically thin disks would be nearly isothermal, in 
which convection does not contribute to the energy trans- 
port and the radiative diffusion approximation breaks 
down. Here, we adopt a one-zone model for optically thin 
disks. Since the radiation pressure is negligible in the 
optically thin medium, the equation of hydrostatic equi- 
librium is given by 



P 



(10) 



where Pgas = pkT/ ^rup is the gas pressure. The viscous 
heating rate is described as 



VIS ^ J- rin^ ^, 



(11) 



where T,.^ is the shear stress tensor. Using the alpha- 
prescription, Trip = —2apga,sF[ , we rewrite the equation 
(11) as 
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where S = 2pH is the column density. The viscosity pa- 
rameter, a, is given by a = min[1.0,10^(iJ/r)^-^], the same 
as in the method for the optically thick disks. The radia- 
tive cooling rate is written as 



Qr 



ad 



2(tT'^i 



(13) 

where r = kE/2 is the optical thickness of the disk. 
Combining the energy balance equation (Qvis = Qrad) and 
equation (10), we obtain the physical quantities of the 
disks. 

To calculate the mass-accretion rate, we use the con- 
tinuity equation {M ~ const, in r-direction) and the 
angular-momentum conservation low. By assuming a 
torque-free boundary condition and angular momentum 
of the gas to be much smaller at the inner boundary than 
at the outer region, we have 



2lT 



(14) 



Thus, using equation (11) and (14). the mass- accretion 
rate is give by 

Throughout the present study, we set the black hole 
mass to be M = 10"^ Mq, unless stated otherwise. 

3. Results 

In figure 1, we present the sequence of equilibrium 
curves on the M — E (mass accretion rate - surface den- 
sity) plane for r~ lO'^rs with Z/Zq = 0, 10~^, and 1, where 
rs = 2GM /c^ denotes the Schwarzschild radius. Here, the 
thick and thin solid lines indicate the resulting equilibrium 
loci for optically thick accretion disks with r > 3, in which 
the energy is transported via convection as well as radia- 
tive diffusion (see §2.1), and the optically thin isothermal 
disks with r < 0.3 (see §2.2). Our numerical method for 
optically thick disks gives optically thin branches for POP 
III disks when A/<lO~^M0/yr. Although they are un- 
physical branches, we represent them by the dotted lines. 

We find that the optically-thick POP 111 disks tran- 
sit to the optically-thin isothermal disks at around M ~ 
lO~^Af0yr~-'^. In contrast, the POP I disks stay optically 
thick (r > 3), even if the mass-accretion rate is less than 
lO^^Afoyr"^. We find that all curves have two major 
turning points at M ^ lO~^M0yr~^ (lower-right turning 
point) and lO^^A/Qyr^^ (upper-left turning point), where 
the gradients of the curves, dM/dJ^, change their signs. 
In other words, the S-shaped equilibrium curves appear 
irrespective of the metallicity. It is known that the mid- 
dle branch {dAI/d'S < : between the upper-left and 
lower-right turning points) is thermally and secularly un- 
stable, whereas the disk is stabilised in upper and lower 
branches, along which dM /dY, > being satisfied. 
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Fig. 1. Sequence of equilibrium curves of the optically-thick 
accretion disks with t > 3 (the thick solid lines) and the opti- 
cally-thin isothermal disks with t < 0.3 {the thin solid lines) 
for Z/Zq = 0, 10-", and 1. The dotted lines are equilibrium 
loci for POP III disks with t < 3, which is obtained by the 
method described in §2.1, although this method might break 
down there. The dashed line indicates the critical mass-ac- 
cretion rate, L-^jf?, where Le is the Eddington luminosity. 
The adopted parameters are M = 10^ Mq and r = lO^rs. 



Figure 1 indicates that if the mass-accretion rate falls 
within the range lO^^Moyr"! <M KlO-^^Moyr-^ the 
POP 111 disks exhibit the limit-cycle behaviour, as well as 
the POP 1 disks (Mineshige & Osaki 1983), since no stable 
branch exists. Note that the range of the mass-accretion 
rate, which induces the limit-cycle behavior, increases as 
the radius increases. We discuss this point in §4.1. When 
the disk stays in the lower stable branch, the surface den- 
sity at r = 10'*rg increases with time, because the mass 
input rate (mass supplied from outside per unit time) is 
larger than the mass-output rate (mass ejected inward per 
unit time). As a result, the disk evolves along the lower 
stable branch toward the lower-right turning point on the 
M — S plane (quiescent state). After a while, it suddenly 
jumps to the upper branch at the right end of the stable 
branch (lower- right turning point). Since the mass-input 
rate is smaller than the mass-output rate in the upper 
stable branch, S decreases along the upper stable branch 
toward the upper-left turning point (burst state) until it 
jumps back to the lower branch. Thus, quasi-periodic os- 
cillation of the mass-accretion rate onto the black hole is 
caused by such state transitions of the disk, leading to 
oscillations in luminosity. 

Here, we stress that the S-shaped equilibrium curve of 
the POP III disks, as well as the POP I disks, does not 
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appear if we employ the one-zone approach for optically 
thick disks (see figure 4 in Mayer & Duschl 2005). It 
is crucial to solve the vertical structure of the optically 
thick disks, including the convective effect, in order to 
investigate the bursting phenomenon caused by a dwarf- 
nova type instability. 

Although the limit-cycle oscillations could arise both 
in POP I and POP III disks, the POP III disks exhibit 
stronger outbursts than do the POP I disks. The surface 
density at the lower-right turning point is about an or- 
der of magnitude larger in the POP III disks than in the 
POP I disks, whereas the upper stable branch is indepen- 
dent of the metallicity. Consequently, the luminosity of 
POP III disks is about an order of magnitude larger than 
that of the POP I disks just after a state transition from 
the quiescent to burst states. As shown in figure 1, the 
equilibrium loci for Z /Zq = and 10~^ almost overlap 
each other. It implies that the magnitude of the outburst 
is independent of the metallicity as long as Z / Zq<\Q~'^ 
is satisfied. Here, we note that the S-shaped equilibrium 
curves for Z / Zq>1Q~^ deviate significantly from those for 
the POP III disks, although we do not plot them in this 
figure. Whereas the upper branch does not depend upon 
Z very much, the lower turning point shifts toward the 
left side as the metallicity increases for Z/Zq >10~^. 

Since we have a ~ 1 and T ~ lO^K in the upper branch, 
the timescale of the burst phase is evaluated to be 



tb-2 



M 



\ IO^Mg 



1/2 



1/2 



yr. 



(16) 



This timescalcs is basically the viscous timescale. The 
timescale for the quiescent phase is given by 7rr^E/M. It 
is described as 



to -4.2 X 10^ 



3 X 1012cm 



M 



3 X 104gcm-2 J I 10-6Afoyr-i 



yr. (17) 



The disc evolves along the lower stable branch on this 
timescale, and it attains the lower-right turning point. 
The timescales of the burst and quiescent phases are much 
longer than those of the transition between the burst and 
quiescent states. Hence, about e = tb/^q 0.5% of the un- 
stable POP HI disks would stay in the burst state, which 
are observed as very luminous objects. Here, we remark 
that global calculations of the disks are need for investi- 
gating the time evolution of the disks in detail. 

In figure 1, we find that the mass- accretion rate exceeds 
the critical rate in the burst phase. In this phase, H /r is 
comparable to, or slightly larger than, unity, although the 
disk is assumed to be geometrically thin in our method. 
The radiation pressure-dominated region would appear at 
the inner region of the disk, in which the disk is unstable 
if the viscosity is proportional to the total pressure. Such 



coupling of two instabilities has been suggested by Lin and 
Shields (1986). By supercritical accretion, the radiatively 
driven outflows would form in the vicinity of the black 
hole, and the photon would be trapped in the flow within 
the trapping radius, rtrap ~ {Mc^ /L^)rs (Ohsuga et al. 
2002). We need a global multi-dimensional approach to 
investigate the disk structure in detail (e.g., Ohsuga et 
al. 2005; Ohsuga 2006; 2007). However, it is beyond the 
scope of the present study. 

The irradiation flux is not taken into consideration in 
the present study. Tuchman et al. (1990) has reported 
that convective disks are stabilized by strong irradiation 
with T > lO^K. However, it is not easy to get the irradi- 
ation flux, since it is very sensitive to the geometry and 
the radiative processes. In the quiescent phase, the disk is 
very geometrically thin, so that the irradiation flux would 
be reduced. In the burst phase, the thick disk forms via 
the supercritical accretion. Then, the radiative flux is col- 
limated in the polar direction, suppressing the irradiation 
effect (Ohsuga et al. 2005). The self-occultation by the 
thick disk would also reduce the irradiation flux (Watarai 
et al. 2005). For this issue, we should study the global 
structure of the disks, fully accounting for radiative trans- 
fer in the multi-dimensional space. 

In figure 2, we represent the resulting optical thick- 
ness of the accretion disks, which are obtained by the 
method described in §2.1, as functions of the mass- 
accretion rate. As shown in this figure, the disks are very 
optically thick in the upper stable branch (burst state), 
Af >10~^MQyr~i. Hence, the spectral energy distribu- 
tions (SEDs) are though to consist of the superposition of 
the blackbody spectra with various temperatures at the 
disk surface. We show them in §4. 

In contrast, as mentioned above, the POP HI disks are 
optically thin in the regime of M <lO~^M0yr~i. Such 
results are unphysical, since our method described in §2.1 
is available only if the disks are optically thick. The 
POP I disks are optically thick, even in the case of 
M <lO-'^M0yr-i. Cannizzo and Wheeler (1984) reported 
that POP I disks become optically thin for a small mass- 
accretion rate in the case of a^O.l, and disks with smaller 
a stay optically thick down to a far lower mass-accretion 
rate. Since a is around 0.03 in the present work, our re- 
sults are consistent with their study. 

In the middle unstable branch, where lO~^M0yr~i <M 

<10"^MQyr~i is satisfied, the optical thickness increases 
very rapidly with an increase of the mass-accretion rate. 
Hence, even in the case of POP HI disks, the optical thick- 
ness is very large in the upper stable and middle unstable 
branches. 

It is found that the spiky structure appears around M = 
lO-^-^Meyr-i for Z/Zq = lO'^. The feature is due to 
the line opacity for heavy elements. Such opacity is very 
sensitive to the temperature as well as the density. Since 
the matter does not contain the heavy elements for Z — 0, 
and since dust cooling is dominant over line cooling for 
Z = Zq, the spiky structure appears only in the case of 
intermediate metallicity. 
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Fig. 2. Optical thickness of disks for Z/Zq = 0, 10"'*, and 
1 as functions of the mass-accretion rate, where we use the 
method described in §2.1. Our method, in which we assume 
the disk to be optically thick, breaks down for POP III disks 
with A'/<lO-^M0yr-^ 

4. Discussion 

4.1. SED of POP III disks 

Wc have so far discussed S-shapcd equilibrium curves, 
which predict the potential possibility of limit-cycle be- 
haviour in POP III disks via thermal and secular insta- 
bilities. We are now ready to discuss whether the limit- 
cycle oscillations are activated in realistic environments 
surrounding POP III disks, or not. 

If we consider black holes in halos with Mhaio ~ 10^ Mq 
at z ~ 20, the virial temperatures of these halos are typ- 
ically T ~ 10'*K, from which the gas is not evacuated by 
radiative feedback of the progenitor star (Kitayama et al. 
2004). 

The density of the gas surrounding the black holes 
would be comparable to, or larger than that of the virial- 
ized halos (Kitayama et al. 2004). Thus, in the first place, 
we assess the mass-accretion rate from the host halo to the 
accretion disk by the Bondi accretion rate with the given 
virialized density/temperature of the halo. 

In figure 3, by using the method described in §2.1, we 
plot the resulting sequence of equilibrium curves of disks 
for r = lO^rs, 3 x lO^rs, and lO^'rs with Z/Zq = 10"* {the 
solid lines). Here, the disks are optically thick, r > 3, on 
the thick solid lines. In contrast, our method described in 
§2.1 gives the optically thin solutions below the lower turn- 
ing points (see the thin solid lines), although this method 
is available only for optically thick disks. We also show the 



Fig. 3. Sequence of equilibrium curves of the accretion disks 
for r = lO'^rg {the lower solid lines), r = 3 X lO'^rg (the middle 
solid lines), and 10*rg (the upper solid lines) in the case of 
Z/Zq = 10~*. Here, we use the method described in §2.1, 
although the disks are assumed to be optically thick in this 
method. The thick and thin solid lines are equilibrium loci for 
the disks with t > 3 and r < 3, respectively. Two dotted lines 
indicate the Bondi accretion rates for p(oo) = 10~^'^gcm~^ 
(upper) and 10~^''gcm~^ (lower). The critical accretion rate 
is represented by the dashed line. The adopted other param- 
eters are M = IO^Mq and Cs(oo) = lOkms-^. 

Bondi accretion rates for p{oo) = 10~^*gcm~^ (the lower 
dotted line) and 10~^^gcm~^ (the upper dotted line). We 
remark that p(oo) — 10~^"'gcm~^ is the typical gas den- 
sity of halos formed at 2 ~ 10 — 30 in standard A CDM 
cosmology. Here, we also assume the sound velocity of 
Cs (00) = 10 kms^^, which corresponds to the halo virial 
temperature, lO^K. As shown in this figure, the S-shaped 
curve shifts toward upper right as the radius increases. It 
is obvious that these accretion disks would exhibit a limit- 
cycle behaviour in the regions of 10"^rs — lO^rs, since the 
Bondi accretion rate corresponds to the middle unstable 
branches at such regions. 

We also investigate the possibility to detect the POP III 
disks. For this purpose, we try to obtain the disk SEDs 
in the burst state. Since the disks are very optically thick 
in this state (see figure 2), we calculate the SEDs by the 
superposition of the blackbody spectra with various tem- 
peratures at the disk surface. Here, we consider the limit- 
cycle oscillation induced by disk instability at r ~ 10'*rg. 
We set the mass-accretion rate to be 2.1 x lO~^M0yr~^ for 
r > lO^'^rs, which is a feasible rate if we assume Bondi ac- 
cretion from a small dark halo. The mass-accretion rates 
within 10'*rs are assumed to be M = 4.5 x lO~^M0yr~^. 
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Fig. 4. SEDs of POP III black hole accretion disks with 
M = W^Mq {the three upper lines) and M = 10'^Mq (the three 
lower lines). Each line type corresponds to z = 7 [the solid 
line), 10 {the dotted lines), and 20 {the dashed lines). The 
SED break denoted by thin dotted lines at lO^'^Hz <i' <10i^Hz 
represents the edge of a Gunn-Peterson trough. 

We can estimate the temperature profile of the disk sur- 
face using above mass-accretion rates and equation (7), 
except for the innermost regions. This equation gives 
the temperature profile of T oc r^'^/^. The mass-accretion 
rate exceeds the critical rate within r = lO'^rg (see fig- 
ure 1). Therefore, photon trapping becomes nonnegligi- 
blc in the innermost regions of the disk, r <{M / LE)rs ■ 
In this region, the temperature profile should be modi- 
fied as T cx r~^/^. Such a profile has been reported by a 
self-similar solution of the slim disk model (Wang & Zhou 
1999; Watarai & Fukue 1999). Watarai et al. (2000) has 
also shown a similar profile by solving the global solution. 

Figure 4 represents the SEDs of face-on POP III disks 
at three rcdshifts: z — 7, 10 and 20. The adopted cosmo- 
logical parameters are SIm = 0.3, fl\ = 0.7, and ho.7 = 1. 
We also add a break due to the Gunn-Peterson effect for 
each SED at lO^Hz <i/ <1015Hz. The resuhing SED shifts 
toward lower left as the redshift increases. 

In this figure, we also plot the SEDs of POP III disks 
around black holes with M = lO^M©. We find that 
the equilibrium curve with M = lO^M© at r = 5 x lO^rs 
roughly overlaps with the curve for M = IO^Mq at r = 
lO^rs. Thus, we consider the limit-cycle oscillation occurs 
at r ~ 5 X lO'^rs. The mass-accretion rates at r < 5 x lO'^rs 
and at r > 5 X lO^rs are set to be M = 4.5 x lO~'*M0yr~^ 
and A/= lO^^A/0yr~^, respectively. We find that the flux 



densities in the case of M ~ 10^ Mq are about 20 — 30% of 
those in the case of M IO^Mq at 10^'''Hz KlO^^Hz. 

Those flux densities are quite low for available facilities 
at present; however, 10 hr integration by OWL telescope 
manages to reach the flux for IO^Mq around 1/im, if they 
could achieve the same limiting flux at 1/im as V-band ^ 

4-. 2. Number of Sources 

Here, we evaluate the expected number density of 
POP III BHs. The number density of POP III BHs 
could be evaluated by the star-formation rate of POP 
III stars. Theoretical studies (Sokasian et al. 2004) pre- 
dict that the comoving star-formation density is roughly 
~ IQ-^Mq yr-i Mpc"^ at ^ >20, although there are stiU 
large uncertainties due to various feedback effects (e.g, 
Susa & Umemura 2006). If we assume all of the POP III 
stars turned out to be POP III BHs, the number density 
of POP HI BHs is roughly evaluated as follows: 

SFR 



np3BH =1-0 X 10-" 
/ t 



10-3MQyr 
/ M 



-iMpc 
-1 



-3 



hljMpc 



(18) 



V lOSyr / V lO^M© , 

Here, t denotes the duration of POP III star formation. 

On the other hand, the comoving volume, which cor- 
responds to the observed field of view with 6 x 6 at 
z ^ z + Az. is 



y =1.6 X 10'' ( — 

>(z) 



10' 
1 



11 



Az 

04 

~i 



/ipyMpC'^ 



(19) 



Here, 4'{z) is a function that equals to almost unity for 
z^lO. Finally, we obtain the observed number of sources 
by multiplying npsBH, V and the efficiency e (see last 
paragraph in §3): 

2 



iVobs =8.0 X lO''^ 



(^) 



1 



10' 

l + z 
11 



Az 

04 

-1 



5 X 10- 



(20) 



Although these are optimistic assessments, it seems pos- 
sible to find 10'^ M0 BHs at z = 10 if we have telescopes 
that cover a 10' x 10' field of view with a sufficient limiting 
magnitude. 

As supplementary evidence, we also evaluate the num- 
ber density of POP III BHs based on the local SMBH den- 
sity. Yu and Trcmaine (2002) assessed the density of mas- 
sive black holes in the present universe using a simple cor- 
relation between the black hole mass and the velocity dis- 
persion of the bulge stars (Gebhardt et al. 2000; Ferrarcsc 
& Merritt 2000). If we consider most of the POP III 
BHs are taken in by massive black holes through merging. 



^ According to the project URL(http:/ /www. eso.org/projects/owl/ 
index_2.html), OWL will have limiting magnitude of V ~ 38. 
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present-day massive black hole density gives the upper 
limit of the POP III BH density. In fact, POP III BHs 
are formed at high-density peaks of a random Gaussian 
density fluctuation of the universe, which tend to be lo- 
cated within density fluctuations with larger scales, such 
as galactic scales. Such intermediate-mass black holes in 
a galaxy settle onto the galactic centre, since they lose 
their orbital angular momentum through dynamical fric- 
tion and three-body interaction processes (Ebisuzaki et al. 
2001). 

According to Yu and Tremaine (2002), the local massive 
black hole density is given by 

PMBH - 2.9 X 105/ig.7MoMpc-^ (21) 

Thus, the upper limit of the comoving POP III BHs num- 
ber density at high redshift is 

np3BH = 2.9 X 10^ (^^^^^ /ig ^Mpc-^ (22) 

It is suggestive that this predicted number density is only 
a few-times larger than the value given in equation (18). 

5. Conclusions 

By solving the vertical structure of the optically-thick 
accretion disks, including the convective energy transport, 
and by employing the one-zone model for optically thin 
isothermal disks, we investigate the stability of the disks 
in primordial environments. 

POP III disks {Z<10-'^Zq) are thermally and secularly 
unstable, and exhibit limit-cycle oscillations, like the POP 
I disks (Z ~ Zq). The outbursts of POP III disks are 
stronger than those of POP I disks. The maximal lumi- 
nosity in the burst state is expected to be an order of 
magnitude larger in POP III disks than in POP I disks. 

The expected flux densities of POP III disks 
surrounding black holes of M = 10"^ M0 are ~ 
10~^^erg s^^ cm~^ IIz~^ aroimd l/.jm in the case of z ~ 
10, whereas the expected number of sources is 0(10"^) per 
10' X 10' field of view. Those disks are within reach of 
future facilities, such as OWL. 
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University. This work was supported in part by a special 
postdoctoral researchers program in RIKEN (KO), and 
by Ministry of Education, Culture, Sports, Science, and 
Technology (MEXT) Young Scientists (B) 17740111(KO) 
and 17740110 (HS). 

References 



Ferguson J. W., Alexander D. R., Allard F., Barman T., 

Bodnarik J. G., Hauschildt P. H., Heffner-Wong A., 

Tamanai A., 2005, ApJ, 623, 585 
Hawley J. F., Balbus S. A., Stone J. M., 2001, ApJ, 554, L49 
Heger, A. & Woosly, S.E. 2002, ApJ , 567, 532 
Hoshi, R. 1979, Progress of Theoretical Physics, 61, 1307 
Huang M., Wheeler J. C, 1989, ApJ, 343, 229 
Kato, S., Fukue, J., & Mineshige, S. 1998, Black-Hole 

Accretion Disks (Kyoto: Kyoto Univ. Press) 
Kitayama,T., Yoshida, N., Susa, H. & Umemura, M., 2004, 

ApJ, 613, 631 
Lin, D. N. C, & Shields, G. A. 1986, ApJ, 305, 28 
Machida M., Matsumoto R., Mineshige S., 2001, PASJ, 53, LI 
Mayer M., Duschl W. J., 2005, MNRAS, 356, 1 
Meyer, F., & Meyer-Hofmeister, E. 1981, A&A, 104, LIO 
Meyer F., Meyer-Hofmeister E., 1984, A&A, 132, 143 
Mii H., Totani T., 2005, ApJ, 628, 873 
Mineshige, S., & Shields, G. A. 1990, ApJ, 351, 47 
Mineshige S., Osaki Y., 1983, PASJ, 35, 377 
Mineshige S., Wheeler J. C., 1989, ApJ, 343, 241 
Omukai K., Nishi R., 1998, ApJ, 508, 141 
Ohsuga, K. 2006, ApJ, 640, 923 
Ohsuga, K. 2007, ApJ, 659, 205 

Ohsuga K., Mineshige S., Mori M., Umemura M., 2002, ApJ, 
574, 315 

Ohsuga K., Mori M., Nakamoto T., Mineshige S., 2005, ApJ, 
628, 368 

Paczyhski, B. 1969, Acta Astronomica, 19, 1 

Pringle J. E., 1976, MNRAS, 177, 65 

Rees M. J., 1984, ARA&A, 22, 471 

Ricotti M., Ostriker J. P., 2004, MNRAS, 352, 547 

Sliakura N. I., Sunyaev R. A., 1976, MNRAS, 175, 613 

Shibazaki N., Hoshi R., 1975, PThPh, 54, 706 

Sokasian A., Yoshida N., Abel T., Hernquist L., Springel V., 

2004, MNRAS, 350, 47 
Smak, J. 1982, Acta Astronomica, 32, 213 
Susa H., Umemura M., 2006, ApJ, 645, L93 
Tuchman, Y., Mineshige, S., & Wheeler, J. C. 1990, ApJ, 359, 

164 

Wada K., Venkatesan A., 2003, ApJ, 591, 38 
Wang, J.-M., & Zhou, Y.-Y. 1999, ApJ, 516, 420 
Watarai K.-y, Fukue J., 1999, PASJ, 51, 725 
Watarai, K.-y., Fukue, J., Takeuchi, M., & Mineshige, S. 2000, 
PASJ, 52, 133 

Watarai, K.-Y., Ohsuga, K., Takahashi, R., & Fukue, J. 2005, 

PASJ, 57, 513 
Yu Q., Tremaine S., 2002, MNRAS, 335, 965 



Abel T., Bryan G. L., Norman M. L., 2000, ApJ, 540, 39 

Balbus S. A., 2003, ARA&A, 41, 555 

Cannizzo J. K., Wheeler J. C., 1984, ApJS, 55, 367 

Ebisuzaki T., et al., 2001, ApJ, 562, L19 

Ferrarese L., Merritt D., 2000, ApJ, 539, L9 

Gebhardt K., et al., 2000, ApJ, 539, L13 

Bromm V., Coppi P. S., Larson R. B., 2002, ApJ, 564, 23 



